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ABSTRACT 

Two  similarity  laws  for  a thick,  axisymraetric,  turbulent  boundary 
layer  on  a long  circular  cylinder  have  been  established.  The  validity 
of  the  assumption  of  constant-stress  moment  in  the  law-of-the-wall 
region  is  analyzed  under  the  assumption  of  similarity.  A new  loga- 
rithmic mixing  length,  which  takes  into  consideration  the  effect  of 
transverse  curvature  on  turbulence,  is  proposed  for  the  law-of-the- 
wall  region.  With  this  logarithmic  mixing-length  model,  the  mean-flow 
momentum  equation  yields  a law-of-the-wall  relation  in  terms  of  the 
exponential  integral  in  the  logarithmic  portion  of  the  inner  layer. 
Comparisons  with  available  data  and  other  proposed  relations  reveal 
that  this  logarithmic  law  gives  a definite  improvement.  Also  a mathe- 
matical explanation  is  given  as  to  why  the  cylinder  data  do  not 
deviate  appreciably  from  the  classical  two-dimensional  law  of  the  wall. 


even  though  it  is  observed  that  the  turbulent  characteristics 


such  as 


the  size  of  large  eddies  and  the  turbulent  intensities  are  much  smaller 


than  for  a flat  plate 


When  the  boundary  layer  is  very  thick  relative  to  the  radius  of 


the  cylinder,  the  flow  in  the  outer  region  of  the  boundary  layer  is 


similar  to  an  axisymmetric  wake.  With  this  assumption,  a velocity 


turbu 


lent  wake.  It  is  assumed  that  the  velocity  defect  has  a separable  form 


2 


Depending  on  the  modeling  of  the  turbulent  shear  stress,  the  mean-flow 
momentum  equation  in  conjunction  with  the  similarity  assumption,  gives 
different  differential  equations.  An  eddy-viscosity  model  leads  to  a 
second-order  linear  differential  equation,  solutions  of  which  are 
known  as  Rummer  functions.  On  the  other  hand,  the  mixing-length  model 
gives  a nonlinear  differential  equation  for  which  the  solution  is 
obtained  in  terms  of  semi-convergent  infinite  series.  The  comparisons 
of  the  results  with  the  available  data  indicate  that  the  eddy-viscosity 
model  is  preferable  to  that  of  the  mixing  length.  The  former  is 
applied  to  compute  the  boundary-layer  characteristics,  and  the  results 
are  shown  to  be  in  good  agreement  with  the  available  data. 
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I NTRODUCT I ON 

A turbulent  boundary  layer  on  a transversely  curved  surface  has  a 
different  structure  from  that  of  a flat-plate  boundary  layer  If  the 
boundary-layer  thickness  6 is  of  the  same  or  higher  order  than  that  of  a 
typical  transverse  dimension  of  the  body.  This  transverse-curvature 
effect  is  of  importance  in  computations  of  the  boundary  layer  near  the 
tail  of  a body  of  revolution,  in  estimating  the  surface  resistance  of 
ship  models,  or  other  three-dimensional  elongated  bodies.  Since  thin 
boundary -1 ayer  approximations  are  not  valid,  these  so-called  thick 
boundary  layers  need  special  treatment. 

The  growth  of  the  turbulent  boundary  layer  on  a circular  cylinder 
of  constant  radius  is  a special  case  of  the  problem  of  transverse- 
curvature  effect,  since  the  pressure  gradient  is  almost  zero  every- 
where and  there  is  no  strong  interaction  between  the  boundary  layer 
and  the  external  flow.  This  problem  has  been  studied  in  the  past  both 
by  experimental  and  theoretical  methods,  because  of  many  engineering 
applications  in  hydrodynamics,  aerodynamics,  instrumentation  and 
textile  engineering,  and  also,  because  of  the  simple  geometry  which 
isolates  the  effect  of  transverse  curvature  from  that  of  the  longitu- 
dinal radius  of  curvature  of  other  elongated  axisymmetric  bodies. 

Although  the  effect  of  transverse  curvature  In  laminar,  axi- 
s>,  metric  boundary  layers  has  been  well  established  by  Seban  and  Bond 


. T. 


investigated  the  transverse-curvature  effect  by  using  the  1/7-th  power- 
law  velocity  profile  and  the  Blasius  skin-friction  law.  He  showed  that, 
at  a given  value  of  momentum- thickness  Reynolds  number,  the  skin  fric- 
tion on  the  cylinder  was  greater  than  that  of  a flat  plate,  and  that 
the  boundary-layer  thickness  was  correspondingly  smaller. 

Most  of  the  investigators  in  this  area  were  concerned  only  with 
the  inner  layer,  particularly  with  the  logarithmic  portion  of  the 
law  of  the  wall.  They  all  assumed  that  the  classical  two-dimensional 
law  of  the  wall  had  to  be  modified.  Some  of  them  tried  to  modify  the 
constants  of  the  two-dimensional  logarithmic  law,  and  others  modified 
the  argument  of  the  function, but  not  the  function  itself.  or  example, 
in  1957,  Richmond  [5]  obtained  a law  of  the  wall  for  the  axisymmetric 
boundary  layer  by  using  Coles'  streamline  hypothesis  in  the  region  near 
the  wall  where  he  assumed  that  mean  flow  is  lomlnated  by  the  wall. 

Richmond  proposed  an  expression  for  the  nond  i.mensional  velocity 

* , * 
u = u /u'f  as  a function  of  the  argument  y (1  + y/2a).  He,  then, 

obtained  the  logarithmic  law  simply  by  replacing  the  argument  of  the 

two-dimensional  logarithmic  law  by  this  new  argument.  In  1967,  Rao 

[6]  showed  that  the  viscous-sublayer  relation  is  not  linear  for  axi- 

•k  k £ * 

symmetric  flows,  but  is  given  by  u = a £n  — . With  y replaced  by 

>V  f >>  { >'< 

a 7,n  • in  u = - ?n  y -f-  B,  he  proposed  a different  form  for  the 


r 


3 

law  oE  the  wall,  La  which  ha  kept  the  constants  of  the  two-dimensional 
law.  Later,  however,  the  experiments  of:  Rao  and  Keshavan  [7]  showed 
that  the  universal  constants  of  a flat-plate  boundary  layer  are 
"functions"  of  Re  and  a for  axisynimetric  flows.  On  the  other  hand, 

il 

Yu  [13]  and  Chin  et  al.  [9]  kept  the  form  of  the  two-dimensional  loga- 
rithmic law  but  tried  to  modify  the  constants. 

Other  work  done  in  this  area  employed  closure  models  of  two-dimen- 
sional flows  and  the  assumption  of  constant  stress  moment  in  the  law- 
of-the-ual.l  region.  In  1963,  Sparrow,  Eckert  and  Minkowycz  [10] 
analyzed  the  turbulent  velocity  and  thermal  boundary  layers  on  a 
cylinder  of  constant  radius  by  using  Deissler's  eddy-viscosity  model. 

In  1972,  Patel  [11]  obtained  a solution  for  the  inner  layer  by  using 
a two-dimensional , mixing-length  model  given  by  Landweber  and  Poreh 
[171. 

fable  I. I summarizes  the  velocity  distributions  in  the  logarithmic 

portion  of  the  law  of  the  wall  proposed  by  various  authors. 

Although  the  boundary  layer  on  a long  slender  rod,  when  6/a  >>  1, 

Is  almost  all  a wake -like  flow,  very  little  attention  has  been  paid  to 

the  ’be  portion  of  the  boundary  layer.  Rao  and  Keshavan  [7]  tried 

to  find  a similarity  variable  for  the  velocity  defect  region  by  trial 

U - u 

and  error.  They  found  that  the  vel.oci.Ly  defect yields  simi- 

* (a  + y)u  ur 

larity  in  terms  of  the  variable  r = , based  on  their 

v 

e pi  i ' ratal  data.  On  the  other  hand,  Yu  [8]  and  Chin  et  al. . [9] 
prt  ni.d  their  ;..vv;iirements  in  classical,  two-dimensional  defect-law 
coordinates  and  t.  r I ■ d to  adjust  the  constants. 


Tn  1972,  White  [1.3]  approached  the  problem  with  a simple  integral 
analysis  in  which  ho  assumed  that  the  law  of  the  wall  suggested  by 
Rao  [6]  is  valid  throughout  the  boundary  layer.  Thus,  the  wake  portion 


ot'  the  boundary  layer  was  not  considered. 

In  1970,  OebecL  [14]  solved  the  turbulent  boundary-layer  equations 
by  an  Implicit  fin ite-diEfercnce  method,  applied  alter  the  momentum 
equations  had  been  linearized.  in  his  treatment,  Cebeci,  used  an 
eddy-viscosity  nodel  oE  thin,  two-dimensional  boundary  layers  which 
ignores  the  eEEect  oE  transverse  curvature  on  turbulence.  However, 

Cebeci's  solution  is  the  only  complete  one  available. 

Tn  1976,  AEzal  and  Narasimha  [15]  studied  the  problem  at  large 

values  oE  the  frictional  Reynolds  number  based  on  the  radius  of 
* auT 

cylinder  a,  a = — , with  the  boundary-layer  thickness  6 of  order  a. 

They  used  the  equations  of  the  mean  flow  and  the  method  of  matched 
asymptotic  expansions  to  show  that  the  flow  can  be  described  by  the 
inner-  and  outer-layer  concepts  that  are  used  in  two-dimensional 
turbulent  boundary  layers. 

Some  of  the  researchers,  however,  did  only  experimental  work; 
namely:  Yanahara  [16],  Willmarth  and  Young  [17]  and  Willmarth,  Wtnkcl, 

Bogar  and  Sharma  [18].  But  the  only  measurement  of  turbulent  quanti- 

! I 

ties  (except  the  wall-pressure  fluctuations  measured  by  Willmarth 
ot  at.  [17,13])  are  given  by  AEzal  and  Singh  [19]. 


Althoug  a large  nu  her  of  measurements  for  l he  mean  velocity 
profiles  are  available,  the  similarity  laws  for  the  axisyimnctric  tur- 
bulent  flow  along  a circular  cylinder  have  not  as  yet  been  well 


r jr- 


considor.it tons,  in  Chapter  2,  a similarity  law  has  been  derived  for 
the  velocity-defect  region,  analogous  to  the  axisymmetric  turbulent 
wake,  both  by  using  eddy -viscosity  and  mixing-length  concepts.  The 
results  are  compared  with  the  data  of  Willmarth  ot  al . [18]. 

In  Chapter  3,  using  the  results  of  the  law-of-the-wall  and  the 
velocity  defect  law,  an  approximate  method  has  been  presented  for 
computing  the  boundary-layer  characteristics. 


sublayer  relation 


For  a Thick  boundary  layer  growing  on  a long,  slender  cylinder  of 
constant  radius,  due  to  the  simple  geometry,  it  is  generally  assumed 
that: 


L)  The  mean-flow  streamlines  remain  almost  parallel  to  the  sur- 
face regardless  of  the  relative  boundary-layer  thickness.  This 
assumption  implies  that  the  normal  component  of  the  mean  velocity 
is  much  smaller  than  the  longitudinal  component,  that  is  v <<  u, 
and  hence  the  induced  pressure  gradient  is  negligible, 
ii)  The  inagnit.  des  of  Reynolds  stresses  are  similar  to  those  of  a 
flat-plate  boundary  layer. 

Under  these  assumptions,  the  boundary  layer  develops  in  a constant 
pressure  field.  In  terms  of  the  cylindrical  coordinates  (x,y)  shown 
in  Fig.  1,  where  x and  y are  the  coordinates  parallel  and  perpendicular 
to  the  axis  of  symmetry  of  the  cylinder,  respectively,  the  boundary- 
layer  cqun , ions  for  the  mean  flow  are 


e 
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whore  u : u(x,y)  is  the  mean  velocity  component  in  the  x-direction. 

v - v(x,y)  is  the  mean  velocity  component  in  the  y-direction. 

r = a 1 y,  with  a the  radius  of  the  cylinder, 

x - i( x , y ) is  the  total  shear  stress,  which  is  given  by 


x = p - • - pu'v' 

3y 


(3) 


Here  p and  p are  the  density  and  the  dynamic  viscosity  of  the  fluid, 
and  pu'v'  is  the  so-called  Reynolds  shear  stress  .lue  to  the  fluctuating 
components  u'  and  v'  of  the  instantaneous  velocity. 

The  boundary  eondil Ions  are 


u ( x , 0 ) = 0 , 


1 1 in  u (x , y)  = U 

y 


v(x,0)  = 0 , 


l (x,  0)  = T (x) 
w 


lim  v(x,y)  = 0 

y >oo 


(4) 


(5) 


where  is  the  velocity  at  i’e  edge  of  the  boundary  layer. 

Since  it  is  assumed  that  the  curvatures  of  the  mean-flow  stream- 


lines are  m-gl  iglble,  the  interaction  between  I he  boundary  layer  and 
the  external  flow  way  be  assumed  to  he  negligible  also.  Therefore, 
once  the  external  flow  is  specified,  Eq.  (!)  and  Kq.  (2)  can  be  solved 
by  using  a proper  closure  relationship  for  the  Reynolds  stresses. 

The  simplest  and  the  most  commonly  used  closure  relations  are  given 
by  the  classical  phenomenological  theories,  which  relate  Reynolds 
stresses  to  the  mean  flow  through  the  well-known  concepts  of  mixing 
length  "fl  '-Idy  viscosity.  Using  the  mixing-length  concept,  Patel  [ll] 
gives  the  velm- ity  distribution  in  the  law-of-the-wall  region  as 


* 


(6) 


x = the  shear  stress  at  y = 0, 
w 

j l - mixing  length. 

For  thin  boundary  layers,  there  are  numerous  formulations  for  the 
mixing- length  distribution.  As  discussed  by  Patel [tl],  the  most  recent 
and  appropriate  one,  given  by  Landwober  and  Poreh[121,  is 


* * / 2 *2 
H = Ky  /t.mh(A  y ' ) 


where  k = 0.418  is  the  Von  K.irmdn  constant, 

A = /J/63  is  a constant  given  by  Patel[ll], 

For  a flat-plate  boundary  layer,  r is  a aimed  to  be  unity  and 
1’qs.  (6)  and  (8)  give  Lhe  law  of  the  wall.  For  an  ax? symmetric,  tur- 
bulent boundary  layer,  however,  It  is  well  known  that  the  total  stress, 
c,  does  not  remain  constant  in  the  law-of- the  vail  region,  although, 
under  certain  conditions,  the  stress  moment  (rr)  can  be  assumed  to  be 

constant  and  equal  to  (ax  ). 

w 


1.2  Analysis  of  Shear-Stress  Oistr ibnl ion 
It  is  usually  assumed  that.  In  a turbulent  flow  along  a smooth 
the  velocity  distribution,  in  the  part  of  the  boundary  layer 


nearest  the  wall,  depends  upon  only  the  density,  viscosity  of  the  fluid 
and  the  wall  shear  stress.  Dimensional  analysis,  applied  to  the  vari- 
ables (n,  y,  a,  x , p,  v) , then  gives  the  functional  relation 
w 


u = f (y  ,a  ) 


(9 


Assuming  that  the  functional  relation  given  by  Kq . (9)  holds  for  a 
thick  axisymmetric  boundary  layer  on  a circular  cylinder  of  constant 
radius,  Rqs.  ( L) , (2)  and  (9)  give  an  expression  for  r as  follows: 


u f(y  ,a  ) - 


= — - f -I-  11 
dx  C 


IL  Y.  ;!llr  M 

_ 3 y ' V 3 a* 


v dx 


JJl  Ay.  , 

3f  3a 

* a x 

. * 3x 

3y 

Da 

i du 

T X 

a (y  f). 

x dx 

_ a y 'c 

The  data  of  VJillmarth  et  al.[l8]  show  that  the  term  Dl'/Da  is  almost 
zero,  as  has  been  verified  analytically  in  the  solution  obtained  for 
the  law  of  the  wall,  which  is  presented  in  a siibsequv  i:  section. 
Therefore  the  above  expression  for  3u/3x  can  be  approximated  by 


Du 

ax’ 


(1 


au 

ay 


2 

ti 

x 

v 


af 

/< 

3y 


(1 


From  Eqs.  (2)  and  (10),  one  obtains 


r! 


or,  1 ntegrati ng  by  parts. 


tin 

v r. 
u dx 

T 


y f - *—5 

a + y 


y E dy 


(12) 


Then  Eqs.  (9),  (10),  (II)  and  (L2)  give 


du 

Dll  , Du  T 

u — + v = u 
3x  Dy  t Jx 


* * ' * f 

£(E+yf)-ffy  + -* 


y £Jy 


a f y ; 


iln 

T 

'x  dx 


c2 

t 

f 

• * 

■y  * A 

y f <iy 

(i3)  ‘ 

- 

a 4-  y 

0 

where  E 


3y 


1-  2 1- 

aiso,  with  x - x c : our  , the  right-hand  side  oE  Eq.  (1)  becomes 
w x 


f 


f, 


3. 

i d ur/v  a •-  * I 

- . (n)  = 3,;  [(a  f y )t 

er  3y  a + y Dy 


(LA) 


(a  + y ) c o " ~ 0 


(a  I-  y )f  dy  + C y Cdy  dy 

j0  J0  J 


* .2  . * 


\y  ry  * * 


o a t:  dy  + y I "dy  I-  i:  y fdy 

J0  'o  '0  J0 


* 

ry  * 2 , 
y f dy5 

0 


Since,  at  y - 0,  r = 1 and  1(0, a ) 0,  the  above  equation  yields 


' ^ * { 
a + y 1 


1 - o'  y f(S)  f(0  l ' - C d^l 

H *j»  * 


For  small  values  of7  y and  a , such  that  the  integrals  in  Eq.  (16) 
are  negl  igible,  t *"  can  be  approximated  by 


+ a 


a + y 


On  the  other  hand,  il:  the  inertia  terms  in  Eq.  (I)  are  neglected  and 
the  remaining  equation  is  ;ntegratod  with  respect  toy,  the  same  expres- 
sion Cor  x"*  , given  in  Eq.  (.17),  is  obtained.  Thus  the  integral  term 
in  Eq.  (16), 


* •>  ' ' (y  r f(v) 

c (y  , a ,a  ) : o f(5)  f(C)  T ,t-  5 d£ 

J0  L a J 


represents  the  effect  of  inertia  on  the  distribution  of  the  shear  stress 


x across  the  boundary  layer  at  a given  section.  The  shear-stress 
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d Lstr  Lbut  i on  given  by  Eq.  (1.7)  has  been  used  widely  in  I he  literature. 
It  assumes  that  the  Inertia  terms  in  the  boundary-layer  equation  are 
negligible,  and  hence  that  the  stress  moment,  rx , is  constant  in  the 
law-of-the-wa 1 1 region,  at  a given  section.  the  relative  magnitude  of 

the  ’inertia’  term  in  Eq.  (16)  can  he  analyzed  as  follows.  In  the 

* * * 
blending  zone,  f(y  ) is  less  than  y , but  of  order  0(y  ).  Then 


P 


f dy 


*3 


* 

y * 


dy 


*3 


Therefore,  from  Eq.  (16), 


+ 

x 


/V 

a + y 


r - 0 


(19) 


This  result  can  be  obtained  directly,  without  using  similarity 
assumptions,  from  the  Taylor-series  expansion  of  the  stress  moment,  re, 
about  y = 0, 


, y a (rx) 

rx  = ax  + ,t  

w 2 2 

dy 


y=0 


y3  a3(tt) 

Dy3 


(7.0) 


y 0 


D(rx) 
since  -- — — 
9 y 


y ”0 


= 0,  as  is  seen  from  the  nous  lip  condition  and  Eq.  (1) . 


Differentiation  of  Eq.  (1)  with  respect  to  y gives 


16 


where  each  subscript  iu  l he  above  equation  corresponds  to  either  a first- 
order  or  second  order  partial  differentiation  with  respect  to  the  sub- 


script letter  used.  From  the  nouslip  condition  at  y = 0,  u = v = u =0, 

x 


from  Eq.  (2)  at  y - 0,  we  have  = 0.  Therefore, 


32(rr) 


3y 


= 0 


y o 


(22) 


Differentiation  of  Eq.  (21),  again,  with  respect  to  y,  with  the  nonslip 
condition  at  y = 0,  and  Eq.  (20),  gives 


.3. 


2u  u + u v = 3-(rxi 
y xy  y yy  pa  3 

3y 


(23) 


r-0 


Differentiation  of  Eq.  (2)  with  respect  to  y and  setting  y = 0 gives 


( u + v \ =0 

l yy)y  0 


(24) 


Also 


y =0 


w 

p 


xy 


y-0 


dr 


w 

dx 


yy 


y 0 


(25) 


aerefore,  from  Eqs.  (23),  (24)  and  (25), 


a3(rr) 


3y 


pa 


2 w dx 


(26) 


y=0 


Substituting  i.qs.  (22)  and  (26)  into  (20)  now  gives  the  Tay lor— series 
expansion  of  the  stress  as 


T = - * T 

a + y w 


UT  T 

+ -2-  -_w  3 

,2  dx  y 


or,  in  nondimens ional  Form, 


a + y 


h 1 ' *3  . 

i 1 - 1 ° y + • • • 


which  is  In  agreement  with  Eq.  (19). 

f 

A typical  value  for  a for  a flat-plate  boundary  layer,  from 

Ref.  [12]  is  a = 6 x 1.0  . Therefore,  at  y = 30,  a--  y = 0.054,  and 

hence,  the  error  in  using  Eq.  (17)  for  the  shear  stress  distribution 

* 

would  be  order  of  5.4%.  As  y in  eases,  the  contribution  from  the 
* * ' 

inertia  term  c(y  ,a  ,o  ) would  increase.  Consequently,  the  error  would 
be  larger  in  using  Eq.  (1.7)  for  the  shear-stress  distribution  rather 
than  Eq.  (16). 


1.3  Numerical  Solution  of  the  Law  of  the  Wall. 

Under  the  assumption  of  similarity,  Eq . (6),  together  with  Eqs.  (8) 
and  (16),  def  ines  the  law  of  the  wall  For  the  axlsyrunetric  turbulent 


boundary  layer.  From  Eqs.  (6)  and  (8), 


2 *2  2 *2  t k 

l [1  + 4k  y tanh(A  y )x  ] 2 


a + y 


' P 

- a 

Jn 


HO  + "V  5 c 


For  various  values  of  a and  o , this  pair  of  equations  has  been 
solved  simultaneously  by  using  the  modified  Adams-Bashfo i Lh  meLliod  and 
Simpson's  1/3  rule,  respectively . The  distribution  of  the  shear  stress 

given  by  F.qs . (16)  and  (1.7)  is  shown  in  Fig.  2.  The  1 aw-of-the-wali 

* * 

range  increases  as  a Increases.  For  the  largest  value  of  y , where 

the  law  of  the  wall  is  valid,  estimated  from  Willmarth's  data,  the 
contributions  of  the  inertia  terms  to  the  shear-stress  distribution  for 

I __  C k 

o = 10  and  for  the  various  values  of  a is  shown  in  Table  1.1. 

TABLE  l.l  ERROR  DUE  TO  INERTIA  TERMS 


* 

a 

k 

y 

P * * 

r = a /(a 
o 

* + 
+ y ) r 

- x+(l  -c) 
o 

. t" 

correction  I o r 

(4) 

50 

169 

0.2283 

0.2265 

0.  794 

100 

2 66 

0.2907 

0.2872 

1.203 

500 

7*7*4 

0.5297 

0. 51.96 

1.893 

1000 

544 

0.6477 

0.6333 

2.217 

Also,  knowing 

* * ' 
that  c(y  , a , a 

f 

) is  much 

less  than  t. 

and  writing 

E + 

Eq.  (16)  in  the  form  t = x (1  - c),  then  I rom  F.q.  (6),  one  obtains 


* y 2x  h(l  - d 2)dy 

* m o 

u ~ / 

J0  1 h /t  + 7*  if  T 


+ a 

Co  * , * 

a + y 


This  shows  that,  the  percentage  error  introduced  by  using  t , rather 
than  x+  = r^(L  - c)  , is  reduced  by  a factor  of  two  in  the  nondimensional 


I 
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velocity  distribution,  u . Therefore  the  effect  of  the  inertia  terms 


on  the  shear  stress  and,  consequently,  on  the  velocity  distributions 


' -5 

can  be  neglected  if  a < 10  , which  is  consistent  with  the  assumed 


k k k 

functional  form,  u = f(y  ,a  ).  However,  if  the  inertia  terms  are  not 


negligible,  the  functional  form  for  u should  oe  f(y  ,a  ,a  ), 

i i * 

where  a = a (R  ,n  ).  Therefore,  throughout  the  remainder  of  this 
study,  the  inertia  terms  will  be  neglected  in  the  law-of-the-wal 1 
region.  Hence  the  distribution  of  the  shear  stress  may  be  approxi- 
mated by  Eq.  (17). 


1 .4  Analysis  of  Distribution  of  the  Mixing  Length 
Generally,  the  flow  in  the  so-called  inner  layer  is  subdivided 
into  three  regions,  namely:  a)  viscous  sublayer,  b)  blending  zone, 
and  c)  fully  turbulent  zone.  In  the  viscous  sublayer,  the  Reynolds 
stresses  are  generally  assumed  to  be  negligible  compared  to  viscous 

k V? 

stresses;  hence  the  mixing  length  is  2 = 0.  With  2 =0  and 

k 

+ 


~r  , Eq.  (6)  gives  the  sublayer  relation  as 


a + y 


u = a 2n(l  + y /a  ) 


(30) 


'or  a flat-plate  boundary  layer,  the  viscous  sublayer  relation  is: 


X X 

u = y 


y < 5 


(31) 


Therefore,  Eq.  (30)  represents  a generalization  of  the  linear  relation 
for  nx [symmetric  boundary  layers.  For  the  fully  turbulent  region,  the 
viscous  stresses  are  generally  neglected  in  comparison  with  the 

k k 

Reynolds  stresses,  and  the  mixing  length  is  taken  as  2 = xy  . 
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For  thin  boundary  layers,  there  are  various  mixing-length  models. 

* * / 2 • 2. 

The  mixing-length  model,  £,  - <y  /tanh(X  y“"  ) , given  by  Landweber  and 

Poreh[12],  and  discussed  by  Patelfll],  gives  the  correct  asymptotic 
behavior  and  is  continuous  throughout  tlie  law-of- the  wall  range.  How- 
ever, in  the  case  of  a thick  axisymmetric  turbulent  boundary  layer 
developing  on  a cylinder  of  constant  radius,  it  is  more  likely  that 
there  is  a direct  influence  of  the  transverse  curvature  on  the  turbu- 
lence, and  hence  on  the  mixing  length,  so  that  it  may  not  be  realistic 
to  adopt  mixing-length  models  established  for  thin  boundary  layers,  as 
has  been  done  in  the  literature.  Indeed,  the  experimental  results 
presented  by  Patel,  N'akayama  and  Damian  [20],  for  the  thick  boundary 
layer  near  the  tail  of  a body  of  revolution  showed  that  mixing  length 
is  directly  influenced  by  transverse  curvature,  and  that  the  mixing 
length  decreases  markedly  as  the  boundary-layer  thickness  increases 
relative  to  the  local  radius  of  curvature. 

The  mixing  length  is  generally  interpreted  as  the  average  distance 
over  which  eddies  transport  momentum  without  losing  their  identity. 
Also,  it  is  known  that  the  large  eddies  contribute  mere  to  momentum 
transfer  than  small  eddies.  Hence  the  mixing  length,  may  be  taken 
to  be  proportional  to  the  size  of  the  large  eddies,  particularly  to  the 
dimension  perpendicular  to  the  direction  of  the  main  flow.  Due  to  the 
geometry,  the  circumference  about  the  circular  cylinder  increases 
linearly  with  the  distance  from  the  wall.  Thus,  in  the  inner  layer, 
eddies  on  a circular  cylinder  can  stretch  more  transversely,  and  hence 
are  shorter  in  the  dimension  perpendicular  to  the  direction  of  he  main 


flow  compared  to  the  eddies  of  equivalent  volume  on  a flat  plate. 
Consequently,  it  is  reasonable  to  expect  the  mixing  length  to  be 
reduced  for  turbulent  boundary  layers  growing  on  circular  cylinders. 

In  other  words,  since  the  region  influenced  by  the  shear  force  on 
circular  cylinders  Increases  as  the  distance  from  the  wall,  the 
intensity  of  shear  force  will  be  decreased  compared  to  the  flat  wall, 
and  hence  the  eddy  size  and  consequently  the  mixing  length  will  be 
reduced.  The  experimental  results  of  Singh  and  Afxal[19],  also  confirm 
that  eddy  sizes  about  a circular  cylinder  are  smaller  than  those  of 
a flat  plate. 

k k 

For  a thin  boundary  layer,  the  sublayer  relation  is:  u = y , 

* * 

and  the  mixing  length  In  the  fully  turbulent  region  is:  & = <y  . 

For  an  axisymmetric  turbulent  boundary  layer,  the  sublayer  rel.ati  t 

k k k k 

is:  u = a 7n(l  + y /a  ).  Hence,  this  suggests,  by  direct  analogy, 
that  in  the  fully  turbulent  region,  the  mixing  length  should  be  re- 
placed by 


l = xa  £n(l+y/a) 


This  analogy  is  reasonable,  because: 

i)  For  a thick  boundary  layer,  Eq.  (32)  yields  smaller  values 
of  the  mixing  length  than  for  a thin  boundary  layer,  as 
required. 

Li)  It  has  the  correct  asymptotic  behavior,  since 


9 j m a £n  (1  + y 


V)  - tin  «*|4  - i(4)  1"  y* 

(a  \a  / ) 
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and  hence  ?.im  x a 7.n(l  + y /a  ) = xy 


Mi)  Mm  a"  ?n  ( 1 i Y ...  ) = Mm  in  y )- 

a"  0 \ a / a"  ><»  a 


’ ?,iin  * 
a"  1 + a y 


■T-  ~r  = 0 
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This  is  consistent  with  Lho  relaminar izat i on  process  observed  in  axi-< 
symmetric  turbulent  boundary  layers.  Consequently,  for  the  axisymmetric 
turbulent  boundary  layer,  Eq.  (8)  will  be  assumed  to  be  given  by 


i\  ,x  ;V  f 2 * 2 

Si  = Ka  £n(l  + y /a  ) ytanh(X  y ) 


1.5  Analytical  Solution  in  the  Sully  Turbulent  Region 
Once  the  distribution  of  the  shear  stress  and  the  mixing  length 
are  established,  the  velocity  distribution  can  be  obtained  from  Eq.  (6) 
by  using  a proper  numerical  integration  scheme.  Substitution  of  Eqs. 
(17)  and  (33)  into  Eq.  ( 6)  gives  the  velocity  distribution  as 


>'c  * 

2a  dy 


n * * \ 2 *2  ? * * 2 * 2 

(a  +y  )<l+  1+4k  a In  (1+y  /a  )tanh(X  y ) 


a + y 


Tn  the  fully  turbulent  region,  from  order-of  magnitude  considera- 
tions, Eq.  ( 34 ) can  be  simplified  and  integrated  analytically.  In  the 

2 *2 

fully  turbulent  region, tanh(X  y ) approaches  unity  and  the  term 


2 *2  2 * . * 
4k  "a  in  (1  + y /a  ) 

Eq.  (34)  reduce  i to 


is  much  greal  r than  one.  Hence 


a T y 


IT  ''  B <a  ) 


* <a  £n(l  + y /a  ) 


ft 


* * 
where  B (a  ) ts  the  Integral  of  Eq.  (34)  from  zero  to  y , the  value 

of  y at  which  flow  becomes  fully  turbulent.  The  change  of  variables, 

■/}  = 1 |-  y '/a',  transforms  Eq.  (35)  to  the  well  -known  Logarithmic 

integral  [21] 


l-  3-  B (a*) 

k I Jnz  r 7 

JrA 

O 


Since  z = J\  + y /a  >1,  Eq. 


•',q.  (36)  may  be  expressed  as  the  well- 


known  Exponential  integral  [21], 


u"  = ~ |E.[2n/l  + y '/a*]  - K.  [9.n/l  + yVa"']  J +■  Br(a*) 


or,  alternatively,  as 


k 1 

u = — E . [ Jin 


A + y/a  ] + B 1 (a*) 


where  B'(a  ) = ^ K ^[ln/l  + y^/a  ] L B^(a  ) 


Also  E . (x)  is  given  in  terms  of  the  infinite  series 


co  n 

r x 


E.(x)  = y + fnx  + 1 — 1 

i n • n! 

n=l 


where  y = 0.5772156  •••  !s  the  Euler  constant.  Therefore,  in  terms  of 
the  infinite  series,  the  velocity  distribution  in  the  fully  turbulent 


region  is  given  by 
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inert. uses  and  y /a  decreases,  the  contribution  from  the  infinite  series 
decreases.  The  following  table  shows  the  contribution  of  the  infinite 


se  ties 

for  some  typ 

* 

Leal  values  of  y 

taken  from  Will  mar th1 

s data  [18] 

TABLE  1 

.2  COMPARISON  OF 

RELATIVE 

MAGNITUDES 

OF 

TERMS  IN 

EQ.  (42) 

a 

y ' 

y /a 

lf 

2^ 

3tt 

4§ 

198.2 

227.4 

1.147 

5.020 

0.3820 

0.4219 

0.73  % 

198.2 

323.5 

1.632 

5.2565 

0.4839 

0.5493 

1.13% 

271.4 

326.3 

1.2023 

5.3672 

0. 3948 

0.4374 

0.73% 

1376 

1981 

1.4393 

7.1123 

0.4459 

0.5009 

0.72% 

* * * 

&n[a  ?,n(l  + y /a  )] 

i * * 

— £n(l  + y /a  ) (first  term  of  series) 

s * n 

y [ £n  ( 1-  + y M l 

n=l  2 n • n! 

% error  due  to  neglecting  the  infinite  series  after  the  first  term. 


Equation  (34)  has  been  integrated  for  various  values  of  X and  com- 
pared with  Willmarth's  data  [13].  The  results  show  that  X is  constant 
inti  equal  to  /3/63.  With  this  value  of  X,  Eq.  (34)  gives  the  complete 
•e loo  tty  distribution  in  the  law-of-the-wal 1 region.  Hence, using  the 
results  from  Eq.  (34),  the  value  of  B(a  ) can  be  calculated  directly 
from  [■>[.  (42) . The  variation  of  B(a  ) with  respect  to  1/a  is  l inear, 


I as  shown  in  Fig.  3,  and  can  be  written  as 

' 

I i 
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B(a  ) = 


41.758 


5.45 


(44) 


Also,  as  shown  in  Table  1.2,  the  largest  contribution  of  the 
infinite  series  conies  from  the  first  term.  Neglecting  the  rest  of  the 
infinite  series,  Eq.  (42)  can  be  written  as 


* 1 
u = — Jin 

K 


,-5.45 

a 


(45) 


For  y /a  <1,  the  argument  of  the  logarithm  can  lie  expanded  as 


(46) 

>v 

Therefore  the  argument  of  the  logarithm, Y , behaves  very  nearly  like 

* * * 
y , even  in  the  neighborhood  of  y /a  = 1.  Also 


* 


3u_ 

* 

3a 


•>'<  * 
3f(y  .a  ) 
* 

3a 


41 . 758 

" *2 
a 


These  results  verify  the  assertion  in  Sec.  1.2  that  3f/0a  is  negli- 


gible. 


28 


[ 

I 
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The  dat.'i  of  Willmarth  et  al  [18]  are  compared  with  Eq.  (34)  and 

with  the  solution  of  Patel  [II]  in  Pig.  4 -Fig.  8.  A comparison  of 

Eq.  (34)  and  Eq.  (45)  is  also  given  in  Fig.  4.  As  can  be  seen  from 

this  figure,  the  analytical  approximation  of  Eq.  (34),  i.e.  Eq.  (45), 

is  in  excellent  agreement  in  the  range  of  the  logarithmic  portion  of 

k 

the  law  of  the  wall,  even  for  the  smallest  value  of  a . Also,  Fig.  9 
compares  the  data  of  Willmarth  et  al  [18]  with  Eq.  (45)  for  various 

values  of  a and  6/a.  It  is  important  to  emphasize  that  the  argument 

k k 

Y given  by  Eq.  (46)  behaves  very  nearly  like  y in  the  neighborhood 

of  y /a  = 1.  This  explains  why  the  Preston-tube  technique  with  the 

usual  calibration  curve  gives  good  measurements  of  the  wall  shear 

stress,  and  why  some  of  the  previous  researchers.  Like  Yu  [8]  and 

Chin  et  al.  [9],  whose  data  are  limited  to  6/a  < 2,  were  able  to  keep 

the  form  of  the  classical  two  dimensional  logarithmic  law. 


I 


d 


Figure  7.  Law  of  the  Wall  for  6/a 


1000  10000 


Figure  9.  Logarithmic  Portion  of  Law  of  the  Wall 
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CHAPTER  2 

THE  VELOC  t lY  - DEFECT  LAW 

2.L  General  Approach 

Unrestricted  rotational  flows  which  are  of  finite  lateral  extent 
and  carried  along  by  the  mean  velocity  of  the  external  free  stream  are 
defined  as  wake  flows.  The  retarded  flow  downstream  of  a solid  body 
placed  in  a free  stream  of  uniform  velocity  is  governed  by  a constant 
total  force  applied  to  the  fluid  by  the  solid  body.  It  is  this,  con- 
stant drag  force  that  determines  the  motion  in  the  wake,  when  self- 
presorving  flow  has  been  attained.  Except  for  the  inner  boundary 
condition,  the  equations  of  the  mean  flow  u and  the  boundary  conditions 
for  an  axlsymmetric  wake  are  exactly  the  same  as  the  mean-flow  equations 
of  a boundary  layer  growing  on  a cylinder  of  a constant  radius,  given 
by  Eqs.  (.1),  (2),  (3)  and  (4)  in  the  previous  chapter.  A thick,  axi- 
symmctric  boundary  layer  is  attained  on  an  infinitely  long  cylinder 
when  the  radius  of  the  cylinder  is  very  small  relative  to  the  distance 
x from  the  nose.  If  the  boundary-layer  thickness  relative  to  the 
transverse  radius  of  curvature  is  large,  the  cylinder  may  be  considered 
as  a small  vorticity-  and  turbulence-producing  disturbance.  Conse- 
quently, the  flow  might  be  considered  similar  to  a wake  flow  with  the 
modified  inner  boundary  condition.  The  important  difference  is  that 
the  drag  generating  the  wake  is  a function  of  the  longitudinal  coordi- 
nate in  the  present  problem.  With  these  considerations,  the 
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velocity  defect  law  has  been  derived  in  a manner  similiar  to  that  which 
has  been  used  for  the  axisymmetric  turbulent  wake  under  the  following 
assumptions: 

i)  As  stated  in  Lhe  previous  chapter,  the  curvatures  of  the 
mean-flow  streamlines  are  negligible;  hence  v <<  u. 

ii)  The  velocity  gradient,  3u/3r,  is  small  in  the  region  of  the 
velocity-defect  law. 

iii)  Oseen's  approximation  is  valid. 

Under  these  assumptions,  the  mean-flow  equation  takes  the  follow- 
ing form: 


e ax 


pr  3r 


where  u , = U - u satisfies  the  boundary  condition 
d e 


im  u (x,r)  = 0 


The  solution  for  the  laminar  boundary  layer  for  the  thick  axi- 

symmetric  flow  along  a circular  cylinder  of  constant  radius,  given  by 

Olauert  and  Lighthill  [ 3 ],  suggests  that  the  proper  axial  length  scale 

2 x 

for  correlating  a thick  hi  alary  layer  is  f;  = -----  , where  Re^  is  the 

a 

Reynolds  number  based  on  cylinder  radius  a,  and  the  free-stream 

velocity  U . On  the  other  hand,  the  solutions  for  laminar  and  turbulent 
e 

wakes  give  the  same  mathematical  form  for  the  velocity  profile  if  the 
eddy-viscosity  concept  is  used.  Therefore,  for  the  present  problem,  the 


use  of  the  following  nondimens ional  quantities  seems  to  be  appropriate: 


18 


.2  x 

^ a Re  ’ 
a 


u , /Re 
d a 


R = 


6 + a 


(49) 


h T 


2 ’ 


r = 


Pu, 


With  these  nondiuensional  quantities,  Kq . (47)  becomes 


02  :'Ud 
F 


9C 


1 3 

~ — (rx  ) 


(50) 


r 3r 


w.  re 


2 Re 


3/2 


a 


U 

e 

a = — 

u 

T 


(51) 


If  the  velocity  variation  across  the  boundary  layer  in  the 
velocity-defect  /-one  is  small  compared  with  the  free-stream  velocity, 
then  the  flow  may  he  self  preserving.  This  suggests  that  the  nondimen- 


slonal  velocity  defect  u^  may  then  be  of  the  separable  form 


»d  = X(0  f(n) 


(52) 


whore 


0 = 


6 + a R 


(53) 


The  boundary  condition  at  n = 1 is 


u , = 0 , t+  = 0 

d 


(54) 


Then  one  obtains 


(55) 


3F, 


= xf(n)  + x 


<lf  dq  dR 
do  dR  dg 


Xf(n)  - Mf’ 


whore 


x(£) 


dx 

df 


R = 


dR 


df 


dC 


dn 


(36) 


and  substitution  of  Eqs.  (52)  and  (55)  into  Eq.  (50)  gives 

6T  ("  ' T n£')  ' - i in  <"I+) 


(57) 


Integrating  Eq.  (57)  with  respect  to  q from  q to  l,  and  using 
boundary  conditions  given  by  Eq.  (54)  yields,  after  some  manipulation, 


(2^  / RX  \ 

nf  dq  ) 

( ' RX  ' J 

n ) 

(58) 


Equation  (58)  is  the  basic  equation  for  the  velocity -defect  law. 

In  order  to  solve  it,  the  variation  of  the  shear  stress  needs  to  be 
modeled.  Depending  on  the  modeling  of  the  shear  stress,  Eq.  (58)  gives 
different  differential  equations  for  x(6)  and  f(q),  as  will  be  shown 
later. 


2. 2 Eddy -Visco sity  A pproach 

The  numerical  solution  of  Cebc  i [14]  showed  that  transverse 
curvature  does  not  affect  eddy  viscosity  directly.  Cebeci  used  the 
free-stream  velocity  and  the  two-dimensional  displacement  thickness  as 
the  characteristic  velocity  and  the  length  scale,  respectively,  in  his 
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treatment.  In  general,  for  free  turbulent  flows  such  as  jets  and  wakes, 
eddy  viscosity,  as  a function  of  the  longitudinal  coordinate,  gives 
very  good  results.  Since  it  has  been  stated  that  the  flow  in  the 
velocity-defect  region  might  be  considered  as  a wake  flow,  the  eddy 
viscosity  may  be  assumed  to  be  a function  of  the  longitudinal  coordi- 
nate, ?,  only.  Hence,  the  distribution  of  the  shear  stress  in  terms 
of  eddy  viscosity,  e(0,  is  given  by 


T = pe 


Substitution  of  the  nondimensional  quantities  given  by  Eq.  (49)  and 
the  velocity  defect  given  by  Eqs.  (52)  and  (51)  into  Eq.  (59)  gives 


2 + e d 3 n 

pu  r = -pe  — -r— 

/Re  3n  3r 
a 


U 

e 

■pe  ~rz 

/Re 


X £. 

aR 


e_  JL 

U2  /Re~  aR 


— R2  J £' 

v R 


Substituting  Eq.  (61)  into  Eq.  (58)  and  rearranging  gives 


V'-  j"2£  + (2i*)  0 


The  similarity  assumption  requires  that 


:(0  e 


2 + R-  = 2b  , 


a constant 


! 

L 


The  solution  of  Eq.  (63),  obtained  by  the  method  of  separation  of 
variables,  is 


X(0  = CR 


2(b-l) 


(65) 


where  C is  a constant  of  integration. 

In  order  to  solve  Eq.  (64) , it  is  necessary  to  know  the  function 


e(^).  For  an  axisymmetrlc  turbulent  wake,  e/v  is  proportional  to  the 

-1/3 


longitudinal  coordinate  x with  a power  law  as  x . Since,  It  has 
been  suggested  that  the  flow  in  the  velocity-defect  region  is  analogous 
to  axisymmetric  wake  flow,  e/v  may  be  assumed  to  be 


e (g)  . m-1 

'v  = H 


(66) 


where  m and  A are  constants.  Hence,  after  substituting  Eq.  (66), 

Eq.  (64)  gives  the  differential  equation  for  the  boundary- layer  thick- 
ness R as 


•*  - & 5* 
a 


(67) 


Integration  of  Eq.  (67)  gives 


/ \ rm+l 

R2  l-  K - % ’ j 

2 m + 1 
a 


A2 


= 2 r a 

a \ in  + l 


(68) 


where  K is  a constant  of  integration. 


From  F.q . (68),  we  obtain 


After  substituting  ^ 


, Fq.  (69)  takes  the  form 
Re  ’ 1 

a 


m+1 


Equation  (65)  then  becomes 


X(C)  = C 


A 


m+1 


2 (b-1) 


or 


m+1 


X = 


(b-1) 


1 + 


(1-b) 


where  B - CA 


2 (b-1) 


From  Eq.  (62),  the  equation  for  f(n)  is 


n f + 2b 


rj  f (1  x]  + a n £ 1 - 0 


Differentiation  of  Eq.  (73)  with  respect  to  q gives 

a2  nf"  + (a2  + n2)C'  + 2(1  - b)  nf  = 0 


lim  f(0  = 0 (76) 

£ ^ 

2.3  So lut ion  o f the  Di f feren t in  1.  Equa t ion  jf  or  f ( nj^ 

The  assumptions  of  similarity  and  the  eddy-viscosity  model  have  led 
to  the  second-order  linear  differential  equation 


£"  (O  + (z  + ^ £’(0  + 2(1  - b)f(c)  = 0 (77) 

if 

where  C = n/a,  and  a and  b are  parameters  to  be  evaluated  from  the 
boundary  condition  and  the  experimental  data.  This  differential  equa- 
tion has  two  singular  points,  at  r,  = 0 and  £ = °°.  Since,  + ~ ) 

2 

= C,  +1  is  analytic  and  r.  gular  • verywhcre,  t;  = 0 is  a regular  sin- 
gular point,  'therefore,  F,q.  (77)  has  at  least  one  nontrivial  solution 
of  the  form 


I 


f(c) 


s 

C 


CO 


k 


6 


(78) 


In  order  to  determine  the  behavior  of  Eq.  (77)  at  infinity,  put 
t = 1/q;  then  the  behavior  of  F,q.  (77)  at  infinity  can  be  obtained  as  t 
goes  to  zero.  With  the  transformation  t = l/s,  Eq.  (77)  then  becomes 


I 


» 


which  s not  regular  at  t = 0.  Since  t ( ^ ) = 1 - V • Therefore, 

V*  tV  t 

t = 0 ai  I hence  5 = » Is  an  irregular  singular  point.  Hence,  a non- 
trivial regular  solution  may  or  may  not  twist  in  the  no ighborhood  of 


C = ■*>  . 

On  the  other  hand,  the  exact  solutions  of  F.q.  (//)  can  easily  be 
obtained  for  the  values  of  b = 0,  1/2,  1. 


Case  1:  b = 0 


f"  + 


+ 2f  = 0 


(80) 


,'hich  can  be  put  into  the  form 


fr  [cf  + r}t)  = 0 


(81) 


Integration  of  Cq.  (81)  once  gives 


Cf*  + 5f  = C. 


(82) 


The  solution  of  this  first-order  linear  differential  equation  is 


-r2/2 

f(C)  = e ^ Z [ C 


1 


,C2/2 


1 d^+C2, 


(83) 


tills  solution  with  = 0 corre  onds  to  the  wake-flow  solution. 


I 
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Case  2:  b = 1/2 


f"  + U + f ' = 


The  solution  of  Eq.  (84)  Is 


2„  / / 2 \ / 2 N 

no  - L i0(v)+c2  K„(v/ 


2 2 
r r 

where  !„(■—)  and  K ( * ) are  the  Bessel  function  of  the  first  kind  and 
4 ' ov  4 

the  modified  Bessel  function  of  the  second  kind,  of  order  zero,  respec- 
tively. 


Case  3:  b = 1 


f"  + (c  + f ' = 0 


Separation  of  variables  gives  the  solution  as 


= C1  . 


4 d/;  + C 


or,  in  scries  form. 


2 4 

£-  + — S_ 

2 2 • i 


This  case  is  impossible,  according  to  Eq.  (65),  since  the  velocity 

defect  u,  decreases  and  R(^)  increases  with  increasing  values  of  £; 
d 

i.e.  b must  be  less  than  1. 


» 
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The  exact  solutions  of  the  differential  equation  for  b = 0,  1/2,  1 
suggest  the  form  of  f(c)  to  be 

HO  = F(z)  , z = £2/2  (89) 

Substitution  of  the  transformation  given  by  Eq,  (89)  Into  Eq.  (77) 
reduces  the  latter  equation  to 

zF" (z)  + (1  - z)  F’(z)  - bF(z)  = 0 (90) 

Eq  ition  (90)  is  a confluent  hypergeometric  differential  equation, 
particularly  known  as  the  Rummer  differential  equation.  It  has  a 
regular  singularity  at  z = 0,  and  an  irregular  singularity  at  z = °°. 

The  method  of  Frobenius  gives  only  one  linearly  independent  solution. 

In  order  to  give  a general  solution  of  the  differential  equation,  it  is 
necessary  to  find  another  linearly  independent  solution.  The  two  lin- 
early independent  solutions  are  given  by  Abramowitz  and  Stegun  [21]  as 

Fi(z)  = 1F1(b,l,z)  = M(b , 1 , z) 

1 ( “ (b)  z” 

F,  (z)  = U(b,  l,z)  = - 777-r  < ,F  (b,l,z)  Inz  + £ — ~~  [<Jj(b +n) 

Z 1 [L  L n=0  (n!) 

- 2*(1  + n)  ] > (91) 

where  ) 

(b)  = b(b  + l)(b  + 2)  ...  (b  + n - 1)  , (b)  = 1 (92) 

n o 

ij/(b  + n)  is  the  logarithmic  derivative  of  the  gamma  function,  which  is 
known  as  the  Digamma  (psi)  function,  and  given  as 


*(*  ' B>  : (n  I)  + b '■  (n  -2)  V b ' 2 ! b + 1 + b + b ''  *<b)  (ri) 


where  i[i(b)  = [j.n  l’(b)]  Is  tabulated  in  Abramowitz  and  Stcgun  [21]. 

F^(b,1,z)  is  the  confluent  hypergeometric  function,  particularly 
known  as  the  Kummer  function.  The  other  notations  for  Kunimer  functions 
are  M(b,l,z)  and  $(b,l,z).  In  series  form,  ^F  (b,l,z)  is  given  by 


. i ^(b, 1 ,z)  £ 2 

L L n=0  (n!) 


in  1 r f (b  + n)  n 

« Z / 1 \ / o •* 


1 y l r i 

= r<b>  nZ0  (n!)2 


Tn  the  neighborhood  of  z = °°,  the  solutions  F^(z)  = M(b,l,z)  and 
F^(z)  = U(b,l,z)  have  simple  asymptotic  expansions  given  by 


“ (1  -b)  (1  - b) 

, / \ ,,/l  . \ b-1  z r 11  11  _“n 

F (z)  = M(b,  l,z)  « z e l j z 

L n=0  11  ‘ 


-l  OT  (b)  (b) 

F (z)  = U(b, 1 , z)  z l ' - (-z) 

£ n=0 


For  z >•  « 


Fx(z)  = M(b,  1 , z)  « -(Lby  e*  zb~L  [1  + 0(]  zl"1)  ] 


F2(z)  = U (b , 1 , z)  * Z_b[l  + Od.r1)] 


Therefore, the  general  solution  for  f(n)  is 

f (n)  = e Z F(z)  = e~Z[C  M(b,l,z)  + C2  U(b,l,z)] 


and,  as  z ->«,  (n  >a>) 


I 


. b-1 

fCn)  * CL  -~y  + C2  e~z  z'b  + 0(|zp) 


(100) 


la  the  outer  region  of  the  boundary  layer,  it  is  assumed  that  the 
mean  flow  behaves  like  an  axisymmetric  wake.  The  solution  of  an  axi- 
symmetric  wake  with  the  eddy  viscosity  assumption  gives  the  defect- 
velocity  profiles  in  terms  of  an  exponential  function.  Hence,  if  the 
boundary  layer  behaves  as  an  axisymmetric  wake  in  the  velocity-defect 
region,  the  constant  must  be  zero.  Therefore,  the  solution  for 

f(n)  is 


f(n)  = e ' U(b, l,z) 


z = -0 — 


( 101) 


From  Eq.  (52),  (65)  and  (101),  the  velocity  defect  is 

nd  = C a2(b_L)  e~Z  U(b, 1, z)  (102) 


where  b,  C and  a are  the  constants  to  be  determined  from  the  experi- 
mental data.  On  the  other  hand,  the  following  analysis  indicates  the 
possible  range  of  the  values  of  b. 

Tiie  nondimens ional  displacement  thickness  for  an  axisymmetric  flow 


is  defined  as 


-r-  d = 


r dr  + 


d - .- 
— r dr 


(103) 


where  r^  is  the  point  of  intersection  of  the  two  laws,  the  law  of  the 
wall  and  the  velocity-defect  law.  Since  the  flow  area  increases  with 
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radial  distance,  the  major  contribution  to  displacement  thickness  comes 
from  the  velocity-defect  region  for  a thick  axisymmetr ic  turbulent  flow. 
Therefore,  the  nondimensional  displacement  thickness  can  roughly  be 
approximated  by 


61  _ fR  Ud  - 

~a  ~ U~  r dr 
Jr  e 
c 


Substituting  Eq.  (102)  and  z 


(104) 


— x— x-  into  Eq.  (104)  gives 
2R  ot 


1 ~ th  a 

— * a CR  e Z U (b , 1 , z) dz  , 

a 

Jz 

o 


Zo  „_2  2 
2R  a 


(105) 


Since  the  displacement  thickness  6^  increases  with  the  boundary  layer 
thickness  6,  b must  be  a positive  number.  Thus  we  have  0 < b < 1. 

From  Eqs.  (91)  and  (94),  the  Rummer  function  U(b,l,z)  can  be 
written  as 


, ( » (b)  z ) 

U ( b , 1 , z ) = -T  T < l [Inz  + i|»(b  + n)-  2<fi(l  + n)  ] > (106) 

lU,;  l n=0  (n!)  ) 


This  series  is  convergent  for  all  values  of  z.  Let  N >>  1;  then 


N-l  (b)  z 

U(b,l,z)  = - l [Anz  + <»(b  + n)  - 2 *(1  + n)]  +Rj, 

n=0  (n!)  (IQ7) 


where 


. . l . y (b>n  * 
% r(b)  <„»* 


[J,nz  + ip  (b  + n)  - 2t|j(l  + n)  ] 


(103) 


'30 


r)L 


Therefore,  Lhe  series  in  Kq.  (106)  can  be  truncated  at  n - 25  for  z < 5. 
For  larger  values  of  z , Llie  asymptotic  expansion  given  by  Eq.  (96) 
should  be  used.  Also,  as  |z|  >■  0,  one  obtains 


U(b, l,z) 


(inz  1 ip (b)  I-  0(|z£nz|) 


(Hi) 


Equations  (70)  and  (101)  contain  six  undetermined  constants; 
namely  A,  m,  K,  b,  C and  a.  These  constants  have  been  determined  from 
the  data  of  Wlllmarth  et  al  [18],  and  Keshavan  [7].  In  Eq.  (69),  as 
£ goes  to  zero,  R goes  asymptotically  to  one.  Hence  the  constant  K 
must  be  of  order  unity.  Consequently,  for  large  values  of  R,  the 


(■  • „» ) 


/ K \ ^ 2 X / 3 

factor  II  + — r J In  Eq.  (70)  is  negligible.  The  plot  of  ( = rs— 

\ R^  / a 

^ -J-  Q 

versus  R = — — on  log-log  paper,  for  large  values  of  R,  then  gave 

cl 

the  value  of  m as  7/13,  from  the  slope  of  the  best  straight  line  drawn 
through  the  data,  as  shown  In  Fig.  10.  Once  the  value  of  in  had  been 
obtained,  the  values  of  K and  A were  evaluated  according  to  Eq.  (68), 
by  using  all  the  available  data.  The  values  found  were  A = 18.5  and 
K = -2. 

The  values  of  b,  a and  C (or  B)  were  determined  solely  from  the 
data  of  Willmarth  et  al  [13].  Let 


S = S (z)  = Ce  U(b,l,z) 

Y • =d  ■ - it)  ^ 


then,  according  to  Eq.  (102), 


Y = S (z)  R 


2 (b— 1) 


(112) 
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Therefore,  for  each  assigned  value  of  /,  = z^,  the  graph  of  Y versus  R on 
log-log  paper  should  give  a straight  line  of  slope  2(b  - 1)  and 
intercept  S(z^),  according  to  Eq.  (112).  Willmarth's  data  were 
plotted  according  to  Eq.  (112)  for  various  constant  values  of  z(or  ri)  , 
as  shown  in  Fig.  11.  The  value  of  b,  determined  from  the  slope,  was 
b = 0.35.  Therefore,  from  Eq.  (72),  x(0  becomes 


■ B>/J  V 


Re  1 - ~ 

a V R2 


(113) 


Hence,  from  Eqs.  (52)  and  (101),  the  velocity  defect  becomes 


U - u i — / x0.65 

-r~-  ■ 8 (>■  - 0 "O’.i-.o  , 


* - 


(114) 


U - u 
e 


-0.65 


(115) 


T = e U(b, l,z)  , 


b = 0.35  , 


2 

z = n~2 
2a 


016) 


After  applying  Eqs.  (115)  and  (116),  Eq.  (114)  becomes 


Z = BT 


(117) 


The  values  of  Z were  obtained  from  Willmarth's  data.  The  values  of  T, 
corresponding  to  the  same  values  of  n , were  calculated  from  Eq.  (116) 
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for  a range  of  values  of  a.  The  value  of  a which  gave  the  best 
straight  line  passing  through  the  origin  was  chosen.  Consequently,  the 
slope  of  the  best  straight  line  gave  the  value  of  B.  The  results  are 
as  follows: 


fri 


18.5 


x/a 


Re 


5/13 


(118) 


U - u 
e 


V-  - Vfe  (‘  - 1) 


0.65 


e-7'  U (b , 1 , /.) 


2a 


(119) 


a = 1/2  , 


b = 0.35 


Comparisons  with  the  data  of  Willmarth  et  al.  [18],  and  Keshavan  [7 
.are  shown  in  Figs.  12  and  13. 


2 . 4 _Tlie_ Mixing  - Le n g t h Approach 

The  mixing  length  for  the  outer  part  of  the  boundary  layer  and  the 
wake  is  generally  assumed  to  be  constant  ana  proportional  to  a charac- 
teristic length  at  a given  section.  For  this  problem,  the  mixing 


Q W i 1 1 ma r t h e t a 1 [18] 


Figure  12. 


Boundary-Layer  Thickness  ( Kddy-Viscosity  Mode  1 ) 


la  terms  of  the  nondimens ional  quantities  given  in  Eqs.  (49)  and 
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(53),  Eq.  (L2i)  becomes 


v U 

+ . 2 e 

T = A ~ 

*•  u Re 
t a 


- \ 2 


(122) 


or,  by  F.qs.  (51)  and  (52), 


t+  = 2A2  /ST  62  X2  f'2 
l a 


(123) 


Equations  (58)  and  (123)  define  the  velocity -defect  law.  Sub- 
stituting Eq.  (123)  into  Eq.  (58)  and  rearranging,  one  obtains 


2A2  /Re 
l a 


q2f  + (l  + 

V RX  / 


qf  dq 


nf 


.'2 


(124) 


The  similarity  assumption  requires  that 


2 A2  /Re~ 
i a 


= Y 


a constant 


(125) 


and 


~ , RX 

2 -1 = n , 

RX 


a constant 


(126) 


Substituting  the  solution  of  Eq.  (126),  which  was  given  previously  by 
Eq.  (65),  with  n = 2b,  into  Eq.  (125)  and  separating  the  variables. 


one  obtains 
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2A2C 

* * “-f*  V^a  ^ 


(127) 


Integration  of  Eq.  (127),  then  gives 


X G ( 3 - n)  „ 

R3_n  10  = — 1 1 jRe  £2 

Y V a 


(128) 


where  (J  is  a constant  of  integration.  Combining  the  group  of  the 
constants  into  single  constant  and  substituting  the  dimensional  quanti- 
ties, Eq.  (128)  yields  the  boundary-layer  thickness  in  the  form 


l 3-n 


R 1 + 


From  Eq.  (65) 


.3-n  i 


- A 


_x/a 


(U 


1 

3-n 


(129) 


= B 


x/a 


2- n 

3- n 


2-n 

v3-n 


1 + 


3-n  , 


( 1 30) 


Also,  from  Eq.  (124),  one  obtains  the  equation  for  f(q), 


f 2 = Ynf  + — 
n 


Of  dq 


(131) 


The  transformation,  f (ri)  =Yf(rl)  reduces  Eq.  (131)  to  one  with  a single 
parameter 


'2  n 

F = qF  + - 
0 


qF  dq 


032) 
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The  boundary  conditions  are 


at  n = 1 , F(l)  = F (1)  = 0 


(133) 


Exact  solutions  of  Eq.  (132)  can  easily  be  obtained  for  the  values 
n = 0 and  n = 2. 


Case  1:  n = 0 


For  n = 0,  Eq.  (132)  reduces  to 


F j— 
P ■ /rl 
/F 


(134) 


Integration  of  Eq.  (134)  with  the  boundary  condition  of  Eq.  (133)  gives 
the  solution 


„ , , 1 ,,  3/2.2 

*0(o)  = 9 <L  “ n ) 


(135) 


This  solution  with  Eqs.  (129)  and  (130)  corresponds  to  that  for 
an  axisymmetric  circular  wake. 


Case  2:  n - 2 


For  n = 2,  Eq.  (132)  becomes 


' 2 2 

OF  = n F + 2 


(1 

nF 

n 


(136) 


Integrating  by  parts  and  using  the  boundary  condition,  Eq.  (136)  yields 


fl  2 u 
n 1 

' n 


(137) 


T .y 


6L 


Dif  feront  i.»  t ion  of  Fq.  (1-37)  gives 


■ ii  i 


F (2qF  + F - p ) = 0 


(138) 


For  F ^ 0,  2qF  + F - q =0  is  a linear  second-order  differential 
equation, of  which  the  general  solution  is 


F2(,l)  = C1  + C2  n'72  + 13  ^ 


(139) 


( 

For  F =0,  F^(n)  = = 0,  according  to  the  boundary  condition  (133), 
an  impossible  solution.  With  the  boundary  conditions  (133),  the  solu- 
tion (139)  becomes 


F2(n)  = — (5  - 6n1/2  + n3) 


I,et  q = (1  - z)  , then 


, . 2 / „ 5 2 2 3 

F (z)  = z ( 1 - 2z  + — z — — z + 
o V 3 3 


1 M 

9 Z ) 


F2(z)  = z 


(!  - f , 


2 2 3 1 

+ z - 3 z + T5 


1 

-5  * ) 


(140) 


(141) 


(142) 


(143) 


Here,  the  factor  z assures  that  the  end  conditions  F(0)  = F (0)  = 0 
are  satisfied. 

The  results  for  n = 0,2  suggest  a solution  for  arbitrary  values 
of  n of  the  form 


F(z)  = l a 


6 2 


Substitution  of  Eq.  (141)  into  Kq.  (132)  gives 


to  a ° i 

F ^(z)  = 4(1  - z)  F - 8n  (1  - z)JF  dz 

•*2 


Differentiating  Eq.  (145),  one  obtains 


--  (f'2)  - 4(1  - z)4  f'  + 8(2  - n)  (1  - z)3  F = 0 
uz 


£ (f'2> 


4 ' V o „P 


(1  - z)  F = l B 


U-z)  F-  l C 

d=2 


Then  Eq.  (146)  yields 


\ = 4 B1 


(145) 


(146) 


(147) 


(148) 


(149) 


(150) 


A = 4 B - 8(2  - n)  C , p = 2,3,4, 

P P P 


(151) 


where,  by  Eqs.  (147),  (148)  and  (149) 


A = T (p  + 1)  m(p  - m + 3)  a a , 
p o m p-m+3 

m=2 


(152) 


'■  ■ \ : 


IS 

p 


(p  + l)a 


P+1 


^P 


a + 6(p  - l)a 
P p-1 


- 4(p  - 2) a 


P-2 


C = a 

P P 


3a  + 3a 
P-1  P-2 


a 

P- 


3 


Then,  from  Eq.  (150) 


and  hence  a^  = 0 or  1.  From  Eq.  (151) 


P+1 

l (p+1)m  (P  -m+3>am  ap-m+3  = 4(P+I)ap+1  “ 8 (2p+2 


+ 24(p  +1  - n)a^  ^ 


8(2p +2  - 3n)a^  2 + 4(p +1  - 2n)a^ 


P = 2,3,4,, 


[f  a2  = 0,  then  from  Eq.  (156),  a^  = a^  = ...  =0  also,  a 
solution.  Hence,  = 1.  From.  Eq.  (156),  for 


P = 2 a3  = -i  (6  - n) 

p = 3 a4  = 3^  (6®  “ 10n  “ a2) 

p = 4 a3  = _SK)  (60  - 16n  + 4n2  - n8) 


(P  -3)a] 


- n)  a 

P 


3 


tr Ivial 


For  p = 6, 7, 8, 9,...,  the  a 's  can  he  successively  obtained  from 
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P+1  4p 

in  - 3 


, P 

— y m(p  - m + 3)a 

n ^ . 


a y (2p  +2  - n)a 

m p-m+3  p(p+l)  p 


+ / ViV  (p  +1  - a) a . - 2 . v-  (2p  + 2 - 3n)a  „ 

p(p+l)  p-l  p(p +1)  p-2 


h - -r-  (p  + 1 - 2n)a 
p ( P PI)  P- 


p-3 


(160) 


Therefore,  the  solution  of  Eq.  (132)  is  the  series  given  by 
Eq.  (144)  with  the  coefficients  given  by  Eqs.  (15/),  (158),  (159)  and 
(160).  For  n = 0 and  n = 2,  this  series  solution  gives  exactly  a 
sixth-degree  polynomial.  For  other  values  of  n in  the  range  0 < n < 2 
the  series  behaves  as  a semi-convergent  one,  which,  if  also  truncated 
at  the  sixth  degree  in  z,  gives  results  in  excellent  agreement  with 
step-by-step  numerical  solutions  of  the  differential  equation. 

Figure  14  compares  the  numerical  solution  of  Eq.  (132)  and  the  series 
solution  with  six  terms  for  various  values  of  n.  Then,  from  Eq.  (160), 
for  p = 5 , 


Thus 

F(z) 


the 


(900  + 168n  - 596n2  + 299n3  - 44n4) 
solution  is  given,  with  sufficient  accuracy,  by 

i - ^ (6  -n)z  + -~r  (60  - lOn  -n2)z2  - (60  - 16n  +4n 

+ (900  + 168n  -596n2  +299n3  -44nSz4 

o L UU  _ 


2 


(161) 


3.  3 
- n ) z 


(162) 


- v *'• 


1 


Mixing-Length  Model  ) 


I 


From  Kqs.  (49),  (52)  and  (130),  wo  obtain 
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Equations  (129)  and  (163)  contain  4 undetermined  constants, 
namely  n,  A,  B and  0.  The  value  of  n = 2b  = 0.7,  was  determined  pre- 
viously in  Sec.  2.3.  Also,  the  value  of  Q_  can  be  estimated  by  using 
the  results  of  the  previous  section.  From  Eq.  (68),  for  small  values 
of  R goes  asymptotically  to  /2,  a value  of  R which  might  be  con- 
sidered to  give  the  virtual  origin.  Therefore,  from  Eq.  (128),  as  £ 
goes  to  zero,  the  value  of  0 is  given  by  0 = -R^  11  = - (/2)^’^  = 2.2. 


The  values  of  Y obtained  from  Willmarth's  data  were  plotted 
against  (lie  values  of  F(z)  obtained  from  Eq.  (162)  with  n = 0.7, 
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corresponding  to  the  same  value  of  z.  The  results  plotted  along  a 
straight  line,  but  the  line  did  not  pass  through  the  origin.  This 
discrepancy  is  due  to  an  ill-defined  boundary-layer  thickness.  Theoret- 
ically, as  n = g -+  ” y 1 (as  y 'y  °°)  , — — - 0.  But,  for  practical 

e _ 

purposes.  Sis  defined  as  the  value  of  y where  -~j — — = 0.01.  There- 

e 

fore,  Eq.  (166)  needs  to  be  corrected  to 
Y = 3 Z + C , 

where  C = constant. 

From  the  slope  and  intercept  of  the  plot  of  Eq.  (166),  the 
values  of  B and  C are 

B = 618 

C = 4.5 

In  summary,  the  results  are 


1 1 
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n = (1  - 7.y 


n = 0.7 


F(z) 


, 53  , 35  2 

1 " 30  Z + 24  Z 


5 3,14 

~9Z  + To  z 


(169) 


Equation  (169)  has  been  obtained  from  Eq.  (162)  by  substituting  n =0.7, 
The  comparison  with  Willmarth's  data  is  shown  in  Figs.  15  and  16, 

The  comparison  of  Fig.  13  and  Fig.  15  reveals  that  eddy-viscosity 
model  yields  a better  correlation  of  the  data  than  the  mixing-length 
model.  Consequently,  the  eddy -viscosity  model  was  selected  in  the 
subsequent  analysis. 


fgure  15  . Boundary-Layer  Thickness  (Mixing-Length  Model) 
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3 . 1.  Boundary-Layer  Character  1st  tcs 
3.1.1  Velocity  Profiles 

Two  similarity  laws  for  a thick,  axisymmetric,  turbulent  boundiiry 
layer  growing  on  a long,  circular  cylinder  of  constant  radius  have  been 
established.  As  for  the  boundary  layer  on  a flat  p1 ate,  these  are  the 
law  of  the  wall  and  the  velocity-defect  law. 

In  Chapter  1,  it  has  been  shown  that  the  logarithmic  mixing 
length,  rather  than  the  linear  one,  gives  a better  modeling  of  the 
Reynolds  tresses  for  axisymmetric  turbulent  flows  with  zero  pressure 
gradient.  Also,  it  has  been  shown  that  the  logarithmic  mixing  length 
gives  the  usual  log  region  in  terms  of  the  Exponential  Integral  of 
£n/IT+  y/a.  This  has  a series  expansion  in  t rms  of  the  argument, 

k k / k >'<  k k k 

Y = a /I  + y /a  £n(l  + y /a  ),  which  behaves  like  y for  y/a  < 1. 

This  result  explains  why  the  cylinder  data,  when  plotted  in  terms  of 
two-dimensional  coordinates,  as  was  done  by  Willmarth  et  al.  [18],  do 
not  deviate  appreciably  from  flat-plate  results.  It  also  verifies  the 
use  of  a Preston  tube  to  measure  the  wall  shear  stress  for  axisymmetric 
boundary  layers. 


As  can  be  seen  from  Fig. 9 , the  logarithmic  region  extends  far- 

k k 

ther  to  larger  values  of  y as  a increases.  It  also  verifies  that 

k 

a = au  /v,  which  is  a measure  of  the  effect  of  transverse  curvature, 

T 
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is  an  important  parameter  for  scaling  the  velocity  distribution  in  the 
law-of-the -wall  region. 

In  Chapter  2,  a si  ilarity  solution  has  been  obtained  analogous 
to  an  axisyuunetr ic  wake  in  the  velocity-defect  region.  The  results 
show  that  the  velocity-defect  profiles  can  be  scaled  by  using  two 
parameters,  x/a  and  R = & + a.  In  order  to  expect  the  boundary  layer 
to  act  like  a wake,  ,$/a  has  to  be  large  compared  to  unity,  which,  in 
turn,  is  a result  of  the  transverse-curvature  effect.  Willmarth's 
data  verify  that  the  portion  of  the  boundary  layer  which  acts  like  a 
wake  becomes  larger  as  6 /a  increases.  The  limiting  case,  that  of  a 
flat-plate  boundary  layer,  (a  ->  «>)  is  thus  excluded  from  consideration. 

In  the  case  of  the  flat-plate  boundary  layer,  the  logarithmic 
parts  of  the  law  of  the  wall  and  the  velocity-defect  law  can  be  ob- 
tained by  assuming  that  there  is  an  over la;  ping  region  in  which  both 
laws  are  valid.  For  the  cylinder,  however,  the  assumptions  made  to 
obtain  the  solutions  fail  at  the  limits  of  the  law  of  the  wall  and 
the  velocity-defect  law.  In  the  law-of-the-wall  region,  it  is 
assumed  that  the  stress  moment  is  constant.  Hence  the  shear  stress 
is  given  by  t+  = /(a+y),  which  ignores  the  inertia  terms  in  the 
mean-flow  momentum  equation.  For  large  values  of  y , however,  the 
inertia  terms  are  no  longer  negligible.  Furthermore,  in  the  deriva- 
tion of  the  velocity-defect  law,  it  is  assumed  that  Oseen's  approxi- 
mation is  valid,  and  that  the  velocity  gradient  3u/Dy  and  the  radial 
velociLy  v are  small.  On  the  contrary,  as  q becomes  small,  both  of 
these  assumptions  fail.  Hence,  the  overlapping  of  the  inner  and  the 


I 


j 
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outer  laws,  which  requires  continuity  of  the  velocity  gradient  Du/3y, 
and  consequently  the  continuity  of  the  shear  stress  which  depends  on 
the  velocity  gradient  Du/3y  through  phenomenological  relations,  does 
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not  occur.  However,  the  velocity  profiles  calculated  from  the  law  of 
the  wall  and  the  velocity-defect  law  do  intersect  each  other  at  an 
angle  which  decreases  as  Re  increases.  Figs.  17-19  compare  the 
velocity  profiles  calculated  from  theory  with  the  Willmarth  data.  As 
can  be  seen  from  these  plots,  the  agreement  is  quite  good  except  near 
the  intersection  point.  Willmarth' s data  reveal  that  the  intersection 
occurs  at  about  y/(6  + a)  = 0.125.  For  all  practical  purposes,  the 
velocity  profiles  can  be  predicted  using  the  law  of  the  wall  and  the 
velocity-defect  law,  and  should  be  smoothed  by  a french  curve  near  the 
intersection  point. 
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3.1.2  Estimation  of  the  Wall  Shear  Stress 

For  the  flat-plate  boundary  layer,  the  overlapping  concept  plays 

an  important  role  in  deriving  a skin-friction  law.  Since  there  is  no 

overlapping  region  for  the  present  problem,  a similar  procedure  for 

obtaining  a skin-friction  formula  is  not  possible.  For  all  practical 

purposes,  however,  the  v ill-shear  stress  can  be  determined  from  the 

intersection  of  the  two  laws,  by  an  iterative  procedure. 

The  velocity-defect  law  gives  the  velocity  profile,  up  to  the 

intersection  point,  by  Eq.  (119).  After  dividing  and  multiplying  by 

a = 0 /u  at  the  intersection  point,  Eq.  (119)  becomes 
ex 

a - u*(r,)  = 4a^  ^1  - j 

I 


0.65  -z 


U(b , 1 , z^) 


(170) 


c© 


n 


where  z = 2r|2  = 2 — ~ 

C C (5  + a)2 


is  the  point  at  which  the  two  laws 


intersect,  indicated  by  the  subscript  c.  From  Eqs.  (37),  (41)  and  (44), 


* Rea 

with  a = , we  have 

a 


Re 

u*(r  ) = - E . [ in/r  ' ] + - in  2 — - + 5.45  - a - 1 

c < x c k a Re  k 


(171) 


Substituting  Eq.  (171)  into  Eq.  (170)  and  rearranging  gives 


oH  I — Zna  = P 

K 


(172) 


where 


H = 1 + - 4 ,a 

Re  \ x 

a 


* (*  ■ #) 


0.65  -z 


U(b,l,z  ) 
c 


(173) 


P = 5.45  + 


— -I  E . [ in/  r ] + in  2Re  - y > 
< l l c a ) 


(174) 


Given  the  values  of  x/a  and  Re  , R = (6  + a)/a  can  be  computed  from 

a 

F.q.  (118)  and  consequently,  z , r , H and  P.  Therefore  a can  be  com- 
puted by  the  iteration  formula 


s (p  - t tn  °i-i) 


(175) 


Table  3.1  compares  the  results  computed  by  using  the  above  pro- 


cedure with  Willmarth’s  measured  values,  and  the  values  of  a obtained 
from  a Clauser  plot. 
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3.1.3  Shear -Stress  Distribution 
The  distribution  of  the  shear 
can  be  approximated  by  t+  = a/ (a  + 
In  the  velocity-defect  region,  the 
of  the  eddy  viscosity  is  given  by 


SL X 


T+  = _ S-  — -A.  f 

v to  W R 


stress,  as  discussed  in  Section  1.2, 
y)  in  the  law-of-the  wall  region, 
shear-stress  distribution  in  terms 
Eq.  (61), 


(176) 


a 


where 


f 


df 

dn 


~ [e  Z U(b , 1, z) ] 


dz 

dn 


e Z U(b,2,z) 


/2z 

a 


(177) 


Substituting  Eqs.  (113),  (118)  and  (177)  and  the  values  of  the  con- 
stants determined  previously  into  Eq.  (176),  one  obtains 


370  i 

13/2  Re 


3/2 


x/a ' 
Re 


29 

26 


1.15 


1 - 


2_ 

„2 


/z 


U (b , 2 , z)  (178) 


Figures  20  and  21  show  the  shear-stress  distribution  computed 
from  Eqs.  (17)  and  (178)  for  two  sets  of  Willinarth's  data.  Around  the 
point  of  intersection  of  two  laws,  there  is  a discontinuity  in  the 
shear-stress  distribution.  As  discussed  in  Section  3.1.1,  this  is 
due  to  the  invalidity  of  the  assumptions  made  in  obtaining  the  solu- 
tions for  the  law  of  the  wall  and  the  velocity  defect  law  in  the 
neighborhood  of  y/(6  + a)  0.125.  For  an  axisymmetr ic , thick  turbu- 
lent boundary  layer,  the  velocity  profiles  are  much  fuller  compared 
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to  those  of  thin  boundary  layers.  The  axisymmetric,  thick,  turbulent 
boundary  layer  profiles  make  an  almost  90  degree  turn  at  about 
y/ (5  + a)  = 0.125.  Thus,  the  classical,  phenomenological  theories  which 
relate  the  shear-stress  distribution  to  the  velocity  gradient  3u/9y, 
through  either  a mixing-length  or  eddy-viscosity  concept,  fail  to  give 
a complete  description  of  the  shear-stress  distribution  if  the  usual 
inner-  and  outer-layer  concepts  are  used.  Apparently,  there  is  a third 
region  between  these  layers  which  would  have  to  be  determined  if  more 
accurate  and  continuous  shear-stress  distributions  are  required. 


3.1.4  Displacement  and  Momentum  Thicknesses 

The  usual  definitions  of  the  displacement  and  momentum  thickness 
for  an  axisymmetric  boundary  layer  flow  are 


and 


0 


1 


a+6 

_u_ 

U 

Ja  e 


(179) 


(180) 


Kelly  [ 2 ],  in  1954,  showed  that  if  the  boundary  layer  is  thick, 

the  displacement  and  the  momentum  thicknesses  require  new  definitions 

\ 

for  cylindrical  flow,  and  the  physically  coherent  definitions  are 


(6  + a) ‘ 


2 

a 


(181) 


(0  + a)  - a - t lL-t)  dr  < 

Ja  e \ e/ 

* 

The  displacement  thickness  6 and  the  momentum  thickness  0 are 
related  to  the  displacement  length  6^  and  the  momentum  length  0^  by 


If  r is  the  intersection  point  where  the  law  of  the  wall  and 
c 

the  velocity-defect  law  intersect,  then  Eqs.  (179)  and  (180)  can  be 
split  into  two  integrals  to  be  evaluated  accordingly, 


r dr  + 


dr 


( 


Also,  Eq.  (180)  gives 


(18! 


(191 


S I nee, 
can  be 


the  term  | 
neglected 


is  very  small  in 


e/ 

in  the  integral  in  Eq. 


the  velocity-defect 
(189).  Hence 


region, 


it 


a 


Tl  + I2 


(19 


Integrating  Eq.  (187)  by  parts  and  using  the  nonslip  condition, 

•k 

u (1)  = 0,  gives 
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1 


1 

2a 


* _ 
u (r  ) 


09 


Si  nee , 

r = 1 +*=  1 +4 
a * 
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* * * 

_3u_  _ 9u  3y 

— k — 

3r  3y  3r 


* 

* 3u 
a — * 
3y 


(19 


Also,  in  the  log-region,  Eq.  (35)  gives 
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Substituting  Eqs.  (193)  and  (194)  into  Eq.  (192),  we  obtain 


rl  * r 

T 11-2*-,.  f 3u  -2  1 f C r ' dr  v,  ,1QRx 

1 20  ‘ C c Jl  3?  K Jin? 


■”  63 

where  r^  = 1 H — , is  the  point  where  the  fully  turbulent  region 
a 

—5/2 

begins.  The  change  of  variables  t = r yields 
*1  = 2a  |?c  U*("c>  - 7 (EiUn?c/2]  " EiUn"l/2])  " 


-l  * 

3u  -2  , , 
— — r dr, 

1 3r 


k _ 

In  Eq.  (196),  u (rc)  Is  given  by  Eq.  (42)  in  Chapter  1,  and  the 

k — 

integral  can  be  computed  by  a suitable  numerical  scheme  with  3u  /3r 
obtained  from  Eq.  (34). 

The  Integrand  of  Eq.  (188)  is  given  by  Eq.  (119).  Hence  sub- 
2 ?2  - - R2 

stituting  z = 2n  =2  , or  r dr=  — - - dz  reduces  Eq.  (188)  to 


■2  _ 2r2 

e z U(b,l,z)dz  , z = — j-  (197) 

z C R 

c 


The  integral  in  Eq.  (197)  is  given  by  Magnus,  Oberhettinger  and  Soni 
[22].  Therefore, 


■ J!  (■  - ?)"  “ ■’  [• 


0.65  « -z  _2 

R e C z ll(b+  l,2,z  ) - 2e  U(b+1 
c c 


,2,2)1 

(198) 


The  integral  1^,  given  by  Eq.  (190)  must  be  evaluated  by  a numeri- 
* 

cal  method  with  u obtained  from  Eq.  (34). 


Table  3.2  compares  the  displacement  and  the  momentum  thicknesses, 
computed  by  using  the  above  relations,  with  Wlllmarth's  data.  In  the 
computations,  Simpson's  1/3  rule  was  used  for  the  numerical  integra- 
tions. The  agreement  is  seen  to  be  good  except  for  the  last  three 
conditions.  A possible  explanation  for  the  discrepancy  for  these 
conditions  is  that  the  flow  was  not  axisymmetric  for  those  cases,  as 
was  indicated  by  Willmarth  et  al.[18]. 
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CHAPTER  4 
CONCLUSIONS 

For  a thick,  ax [symmetric , turbulent  boundary  ’ayer  on  a 
circular  cylinder  forms  of  two  similarity  laws,  one  a law  of 
the  wall,  the  other  a velocity-defeet  law  have  been  derived. 

In  the  region  of  the  law-of-the-wall,  the  logarithmic  mixing 
length  is  a better  approximation  than  the  linear  one  for  the 
Reynolds  stresses  for  thick,  axisymmetric  turbulent  flows.  With 
this  logarithmic  mixing  length,  the  mean-velocity  distribution, 
expressed  as  an  integral,  Eq.  (34),  or  by  an  approximation  to  it 
which  yields  a logarithmic  law  in  terms  of  the  exponential 
integral,  Eqs.  (42)  and  (45)  shows  excellent  agreement  with  the 
data.  Eq.  (45)  is  recommended  as  a simple,  analytical  repre- 
sentation of  the  law  of  the  wall  for  the  thick,  axisymmetric, 
turbulent  boundary  layers. 

* 

The  frictional  Reynolds  number  based  on  cylinder  radius,  a , 

is  an  important  parameter  in  the  law-of-the-wall  region.  For 
* 

large  values  of  a , the  law-of-the-wall  approaches  the  classical 
two-dimensional  form. 

* 

Even  for  the  intermediate  values  of  a , in  the  region  where 
y /a  <1,  the  argument  of  the  logarithmic  law 
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Y = a /l  + y /a  ?,n(l  + y /a  ) behaves  like  y . Therefore,  the 
Pres ton-tube  technique  can  be  used  to  measure  the  wall-shear 
stress  in  a thick  boundary  layer,  with  the  usual  calibration 
curve,  as  was  assumed  by  Willmarth  et  al  [18]. 


2 x/a 

(5)  The  axial  length  scale,  £ = , of  the  thick  axisymmetric, 

a 

laminar  boundary  layer  is  also  a proper  scale  for  correlating  a 

thick,  turbulent  boundary  layer.  With  this  scale,  the  boundary- 

layer  thickness  is  well-correlated  as  power  law  or  a wide  range 

of  values  of  5/a,  x/a  and  Re  . 

a 


(6)  A velocity-defect  law,  in  separable  form  u^  = X(£)  f(n)>  has 

been  obtained.  This  similarity  assumption  yields  the  nondimen- 
sional  velocity  defect,  ^ f(n).  This  result  shows 

that  the  velocity-defect  profiles  depend  strongly  on  5/a,  as 
indicated  by  Willmarth  et  al  [18]. 


(7)  The  eddy-viscosity  concept  for  the  velocity-defect  region 

appears  to  be  preferable  to  that  of  mixing  length,  but  the  form 
of  the  eddy  viscosity  needs  verification  by  turbulent-stress 
measurements.  The  eddy-viscosity  model,  in  conjunction  with 
similarity,  gives  the  velocity  defect  in  terms  of  particular 
confluent  hypergeometric  functions  known  as  Rummer  functions. 

The  inixing-length  model,  on  the  other  hand,  leads  to  a nonlinear 
differential  equation,  the  solution  of  which  is  a semi-conver- 


gent infinite  series. 


■— 
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(8)  The  wake  portion  of  the  boundary  layer  increases  as  6/a 
becomes  larger  and  the  wall  region  decreases  correspondingly. 
Approximately  1/8 -th  of  the  boundary  layer  obeys  the  law  of  the 
wall. 

(9)  The  two  similarity  laws  do  not  overlap.  However,  they 

intersect  at  about  • = 0.125.  Boundary-layer  characteristics 
6 + a 

computed  by  using  this  intersection  point  agree  with  the  data 
within  an  acceptable  range  of  error. 


(10)  There  is  insufficient  data  at  large  values  of  6/a,  for  a 

variety  of  situations,  to  verify  the  accuracy  of  the  velocity- 
defect  law  obtained  in  Chapter  2.  As  stated  by  Willmarth  et  al 
[18], 

"In  the  future  we  may  gain  enough  knowledge  of  the 
effects  of  transverse  curvature  to  consider  the 
possibility  of  creating  an  empirical  formulation 
for  the  mean  flow  in  an  axisymmetric  boundary  layer. 

This  was  only  possible  in  the  two-dimensional  case 
after  sufficient  data  and  understanding  had  accumu- 
lated, ..." 

In  the  present  work,  such  a formulation  has  been  proposed,  but, 
in  view  of  the  paucity  of  the  available  data,  it  may  be  necessary 
in  the  future  to  modify  the  values  of  the  constants  when  addi- 

' i 1 

tional  data  become  available. 


Wtf'  ' • r ' r ' . *»■ 


Cylinder  in  Axial  Incompressible  Flow,"  Journal  of  Aeronautical 
Sciences,  Vol.  21,  1954. 

[3]  Glauert,  M.  B. , and  Lighthill,  M.  J. , "The  Axisymmetric  Boundary 
Layer  on  a Long  Thin  Cylinder,"  Proceedings  of  the  Royal  Society, 
London,  Vol.  230,  1955. 

[4]  Landweber,  L.,  "Effect  of  Transverse  Curvature  on  Frictional 
Resistance,"  David  Taylor  Model  Basin  Report  No.  689,  1949. 

[5]  Richmond,  R. , "Experimental  Investigation  of  Thick  Axially 
Symmetric  Boundary  Layers  on  Cylinder  at  Subsonic  and  Hypersonic 
Speeds,"  Ph.D.  thesis,  California  Institute  of  Technology, 
Pasadena,  Ca. , 1957. 

[6]  Rao,  G.N.V.,  "The  Law  of  the  Wall  in  a Thick  Axisymmetric  Turbu- 
lent Boundary  Layer,"  Journal  of  Basic  Engineering,  Trans.  ASME, 
Series  D,  Vol.  89,  No.  1,  March  1967,  pp.  237-238. 

[7]  Rao,  G.N.V.,  and  Koshavan,  N.  R. , "Axisymmetric  Turbulent  Boun- 
dary Layers  in  Eero  Pressure  Gradient  Flows,"  Journal  of  Applied 
Mechanics,  Trans.  ASME,  Vol.  94,  1972,  p.  25. 

[8]  Yu,  Y.  S.,  "Effects  of  Transverse  Curvature  on  Turbulent  Boundary 
Layer  Characteristics,"  Journal  of  Ship  Research,  Vol.  3,  1958, 
p.  33. 

[9]  Chin,  Y.  T.,  Hulesbos,  J.,  and  Hunnicutt,  G.  H.,  "Effect  of 
Lateral  Curvature  on  the  Characteristics  and  Skin  Friction  on  a 
Turbulent  Air  Boundary  Layer,  With  or  Without  Helium  Addition," 
Proceedings,  Heat  Transfer  and  Fluid  Mechanics  Institute,  1967. 

[10]  Sparrow,  E.  M. , Eckert,  E.R.G.,  and  Minkowycz,  W.  J.,  "Heat 
Transfer  and  Skin  Friction  for  Turbulent  Boundary  Layer  Flow 
Longitudinal  to  a Circular  CyLinder,"  Journal  of  Applied  Mechan- 
ics, Vol.  30,  1963,  p.  37. 


92 


[11]  Patel,  V.  C. , "A  Unified  View  of  Che  Law  of  the  Wall  Using 
Mixing-Length  Theory,"  Iowa  Institute  of  Hydraulic  Research, 
Report  No.  137,  1972. 

[12]  Landwcber,  L. , and  Porch,  M. , "Frictional  Resistance  of  Flat 
Plates  in  Dilute  Polymer  Solutions,"  Journal  of  Ship  Research, 

Vol . 17,  No.  4,  1973,  pp.  227-240. 

[13]  White,  F.  M. , "An  Analysis  of  Axisymmetric  Turbulent  Flow  Past 
a Long  Cylinder,"  Journal  of  Basic  Engineering,  Trans.  ASME, 

Vol.  94,  1972,  p.  200. 

[14]  Cebeci,  T.,  "Laminar  and  Turbulent  Incompressible  Boundary  Layers 
on  Slender  Bodies  of  Revolution  in  Axial  Flow,"  Journal  of  Basic 
Engineering,  Trans.  ASME,  Series  D,  Vol.  92,  p.  545. 

[ L5]  Afzal,  N.,  and  Narasimha,  R. , "Axisymmetric  Turbulent  Boundary 
Layer  Along  a Circular  Cylinder  at  Constant  Pressure,"  Journal 
of  Fluid  Mechanics,  Vol.  74,  1976,  p.  113. 

[15]  Yasuhara,  M.,  "Experiments  of  Axisymmetric  Boundary  Layers  Along 
a Cylinder  in  Incompressible  Flow,"  Transactions  of  the  Japan 
Society  of  Aerospace  Sciences,  Vol.  2,  1959,  p.  33. 

[17]  Willmarth,  W.  W. , and  Yang,  C.  S.,  "Wall-Pressure  Fluctuations 
Beneath  Turbulent  Boundary  Layers  on  a Flat  Plate  and  a Cylinder," 
Journal  of  Fluid  Mechanics,  Vol.  41,  p.  47. 

[18]  Willmarth,  W.  W. , Winke'l,  R.  F,.,  Bogar,  T.  J.  , and  Sharma,  L.  K.  , 
"Axially  Symmetric  Turbulent  Boundary  Layers  on  Cylinders:  Mean 
Velocity  Profiles  and  Wall  Pressure  Fluctuations,"  The  University 
of  Michigan,  Department  of  Aerospace  Engineering  Gas  Dynamics 
Laboratories,  1975. 

[19]  Afzal,  N.,  and  Singh,  K.  P.,  "Measurements  in  an  Axisymmetric 
Turbulent  Boundary  Layer  Along  a Circular  Cylinder,"  Aeronautical 
Quarterly,  Vol.  27,  1976,  p.  217. 

[20]  Patel,  V.  C.,Nakayama,  A.,  and  Damian,  R. , "An  Experimental 
Study  of  the  Thick  Turbulent  Boundary  Layer  Near  the  Tail  of  a 
Body  of  Revolution,"  Iowa  Institute  of  Hydraulic  Research, 

Report  No.  142,  1973. 

[21]  Abramowitz ,M. , and  Stegun.I.A.,  "Handbook  of  Mathematical  Func- 
tions," Dover  Publications,  Inc.,  New  York,  1970. 

[22]  Magnus, W.,  Oberhet t inger , F. ,and  Soni.R.P.,  "Formulas  and  Theorems 
for  the  Special  Functions  of  Mathematical  Physics."  Springer- 
Verlag,  New  York  Inc.,  1966. 


93 


APPENDIX 

+ 

SOME  PROPERTIES  OF  RUMMER  FUNCTION  U(b,c,z) 

The  following  properties  of  the  Kummer  function,  U(b,c,z), 
applied  in  the  text  and  were  required  in  numerical  evaluations. 


,n 

(1)  — --  [e  Z U(b,c,z) ] = (-1)°  e Z U(b,c+n,z) 

dz 


(2) 


e Z zC  ^ U(b,c,z)dz  = -e  Z zC  U(b+l,c+l,z) 


(3)  (b  + z)  U(b , c , z)  + b(c  - b - 1)  U(b  +l,c,z)  -zU(b,c  +l,z)  = 0 


Magnus  et  al. 
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